Influence of four-wave mixing and walk-off on the self-focusing of coupled waves 
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Four-wave mixing and walk-off between two optical beams are investigated for focusing Kerr media. 
It is shown that four-wave mixing reinforces the self-focusing of mutually-trapped waves by lowering 
their power threshold for collapse, only when their phase mismatch is small. On the contrary, walk-off 
inhibits the collapse by detrapping the beams, whose partial centroids experience nonlinear oscillations. 
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Self-trapping of intense electromagnetic waves in non- 
linear media have inspired interest since the earliest days 
of nonlinear optics ||l| . For a focusing Kerr medium, this 
process is often described by a (l-l-l)-dimensional nonlin- 
ear Schrodinger equation (NLS), when, for instance, opti- 
cal pulses are assumed to only undergo anomalous group- 
velocity dispersion (GVD) while their transverse diffrac- 
tion is neglected. Recent developments based on this 
model rapidly increased due to the possibility of making 
two wave components couple nonlinearly and propagate 
as mutually-trapped beams |^|-^. In the seventies, such 
optical beams with two orthogonal polarized components 
were theoretically predicted to form spatial solitons pre- 
serving their incident shapes and powers 1^. However, 
these solitons were only recently detected in experiments 

, partly because of the difficulty in overcoming para- 
metric four-wave mixing (FWM) processes which become 
relevant when the two orthogonal polarizations are not 
incoherently coupled Besides FWM, they moreover 
undergo linear convection, often termed as "walk-ofF', 
which can be either attached to the half-difference of their 
group velocity in birefringent media |^,^, or connected 
with the angle between the transverse and carrier wave 
vectors in, e.g., biased photorefractive crystals 0. 

In the previous works, only one-dimensional dispersion 
was considered and the interaction patterns between two 
partial pulses were analyzed in terms of stable solitons 
keeping a fixed sech-shape ||J^. However, the under- 
lying assumption that the soliton components preserve 
their shape while propagating does not hold at higher 
dimensions. It is indeed well-known that NLS solutions 
are unstable, by spreading out or self-focusing until col- 
lapse, when transverse diffraction is no longer disregarded 
and enters (Z?+l)-dimensional equations with a number 
of spatio-temporal dimensions, D, at least equal to two 

. High-power pulses then sharply increase in amplitude 
and undergo destructive modulations, which significantly 
affect their mutual interactions Even if the collapse 
singularity is ultimately arrested by saturation mecha- 
nisms, the dynamics preceding it is important for under- 
standing this growth stage, which justifies to study the 
self-focusing of nonlinear waves when FWM and walk-off 



cannot be ignored. In this letter, the inffuence of both 
these effects is investigated from (Z)-|-l)-dimensional cou- 
pled NLS equations, with emphasis on the potential mod- 
ifications of the power threshold for collapse when D=2. 
By means of analytical estimates and numerical verifica- 
tions, we show that FWM contributes to self- focusing by 
lowering this threshold for pulses with small phase mis- 
match only, whereas walk-off acts against the collapse by 
inducing oscillations in the wave trajectories. 

We consider the slowly- varying envelopes ui(r, z) and 
U2{r,z) of two waves copropagating along the z-axis of 
a Kerr medium. Here r refers to the vector of a gener- 
alized D-dimensional transverse diffraction plane, which 
can also include a retarded time variable when anoma- 
lous GVD is retained. The waves undergo walk-off and 
FWM as described by the generic model 



+Bu'u;^ = 0, n=l,2, 



(1) 



where * means complex conjugate. Si — —62 — S and 
Pi = —(32 = p. In Eq.(I), standard notations have 
been used. In particular, the second term describes 
walk-off effects, the third term accounts for the mis- 
match in wave numbers between the components and the 
last term represents FWM. The constants A and B are 
positive and they measure the strength of the nonlin- 
ear coupling between the two waves. In isotropic me- 
dia, their values satisfy A -\- B=l when ui and U2 rep- 
resent two orthogonal polarizations of a vector field (as, 
e.g., A=2B=2/3 in birefringent fibers |^), while they are 
linked by A=2B—2 for two copropagating beams with 
scalar amplitudes ui and U2. For Eq.(|^) the total power 
P = Pi+ P2 = \\ui\\l + IIW2II2 is always conserved, and it 
keeps invariant each individual power P„ = ||u„||2 when 
-6 = only. For notational convenience, we make use 
of the standard norms, ||/||^ = / Ifl^dr, where / df 
denotes an integration over the D-dimensional transverse 
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space. Equations (|l|) also conserve the Hamiltonian 

2 

^ = E{ll^^""ll2 + ^nPn - Xn{A,B) 

Tl=l 



where 

Xn{A,B) = 




is the nonhnear potential related to the Kerr and FWM 
contributions. Furthermore, we can derive a virial iden- 
tity describing the evolution of the total mean square 
radius of both waves along z. This relation governs the 
second-order z-derivative of the integral 

[ |r-(r)p^Kpdf, (3) 

n=l 

which involves the total center of mass (f) = 
I ^X]n=i Wnl'^d'i^- After employing the straightfor- 
ward procedure expounded in ||^, the virial identity is 
found to read 

d'^I = dll/P - 2|a,(r)p - 2(r) • dl{r) (4) 
where I ~ j ^Li \un?dr = P{1 + |(r)|2), 

2 

dll = 4DH + 4^{(2 - D)\\V^Un\\l - DPnPn 

n=\ 

+ ^ _ ■ Im 1 r{ul_^u':)dr}, (5) 

2 

Pd, (r) =M' = M" + Y, SnPn, (6) 



M° = 2Im ^ / {u:y^un)dr, (7) 

n=l 

with = Here, denotes the conserved total 

momentum. When collapse develops, the integral I{z) 
vanishes at a finite propagation distance Zc, implying 
In = Pn^\\{^~ (^)^n|l2 — > as z ^ Zc. In this limit, the 
gradient norms of m„ must blow up by virtue of the in- 
equality Pn < (4/Z?^)I„||V_LM„|j2, which makes the wave 
amplitude diverge in turn. In NLS systems, such a blow- 
up generally takes place at propagation distances shorter 
than Zc [||, and it arises when the total power exceeds 
some threshold. For instance, when considering a single 
wave governed by the (2+l)-dimensional NLS equation 



idzU + V\u + = 0, blow-up can only develop if P 
exceeds the threshold value Pc — 11.7, which is close to 
the critical power An fixing H = for Gaussian beams 
|^,|o|. In the present context, we determine the influ- 
ence of FWM and walk-off on self-focusing by dividing 
the analysis into two separate parts devoted to FWM 
alone and to walk-off alone. 

1. FWM without walk-off : For the case B ^ Q and 
5 — the virial relation (^) can be bounded from above 
when 13 > 2 as follows 

dlX < AD{H - f3Pi + f3P2)/P - 2|M°|V^'^ (8) 

where the strict equality applies to the case D = 2 only. 
Thus, a total collapse of the wavefunction, in the sense 
I(z) — > 0, always occurs for D > 2, whenever {H + \(3\P) 
is negative. In (2-f 1) dimensions, the threshold power for 
collapse can be estimated by bounding the Hamiltonian 
from below with a combination of Cauchy-Schwarz and 
Sobolev [Ijunlll < (2/Pc)|| V_LU„||2||u„||2] inequalities : 

H+mP>j2i^-^]\\V^uJl (9) 

n=l c 

Since the left-hand side of Eq.(^ is finite, a wave collapse 
with ||Vj^M„||2 +00 cannot occur when the individual 
power of both components is below the threshold 

= Pe/(l + A + S). (10) 

Hence, P„ < is a sufficient condition for no collapse. 
In terms of the conserved total power P, collapse is surely 
avoided when P is below P^, but due to power trans- 
fers between both components, the bound P < can 
significantly underestimate the actual threshold for col- 
lapse when i? 7^ 0, as was observed in similar studies of 
quadratic nonlinear media . In spite of this, when the 
two components are initially identical, Eq.(P with /3 = 
reduces by symmetry to the single-wave NLS equation 
idzU+V\u+[l+A+B)\u\^u = and power transfers are 
eliminated. thus restores the exact threshold power 
for collapse of each component and, therefore, P > 2V^ 
can here yield a reasonably-good estimate of the exact 
collapse threshold for small mismatch values, provided 
that Mi(f, 0) ~ U2{r,Q). In this situation, Eq.([lO|) shows 
that FMW lowers the power threshold for collapse and 
strengthens the self- focusing dynamics, compared with 
that of incoherently-coupled waves [B — 0) ||] . 

Using an iterative, radially-symmetric, mid-point 
Crank-Nicholson finite-difference scheme, we have per- 
formed numerical simulations of Eq.(|l|) in (2+1) dimen- 
sions with the Gaussian initial condition 

u„(f, 0) = v/^V^exp (-rV2A2). (11) 

We used a resolution of A?- = 10^^ over the interval 
r=[— 20, -1-20] and a stepsizc of Az ~ 10"'^, which kept 
the relative deviation of P from its initial value below 
10^** over a distance z=[0,50]. 
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FIG. 1. Collapse threshold power 27^,? versus B for 
ui(r, 0)=U2(r, 0) with A— 2 (solid line). The dotted line rep- 
resents 2P^. The black and white marked points indicate 
the numerically-calculated collapse thresholds for [3 — and 
/3 = 20, respectively. 



The solid line in Figure 1 shows the collapse threshold 
2V^ versus the FWM coefficient B when Pi = P2 = P/2 
and A = 2. For 13=0, this threshold power agrees with 
its numerical counterpart, since the dynamics is entirely 
determined by symmetry from the single- wave NLS equa- 
tion. For P ^ 0, this symmetry reduction does no longer 
apply and the waves exchange power. However, at large 
mismatches (/? ^ 1), the FWM terms can merely be av- 
eraged out (B 0) from Eq.(|^) after applying the phase 
transform u„ — > u„ exp (— i/?„z), so that power transfers 
disappear and the collapse threshold is simply given by 
27^°, as confirmed by the white marked points of Fig. 
1. In contrast, a more intriguing behavior appears for 
moderate values of /3. For /3 ranking as 1 < /3 < 10, the 
two components were indeed observed to diffract with 
no significant power exchange in the range P < 2V^. 
For P above 2V^, they started to periodically transfer 
their power between each other, which generated strong 
resonances following several focusing/defocusing cycles. 
Rapid oscillations in each partial power created spiky 
resurgences in the wave dynamics, which made it com- 
putationally difficult to determine an accurate collapse 
threshold. Despite this, a non-zero (3 always counter- 
acts the self-focusing promoted at exact phase matching 
{(3 = 0), by increasing the power threshold for collapse. 

2. Walk-off without FWM : We now consider the oppo- 
site situation 5^0 and B = 0. As the nonlinearities only 
depend on intensities, the phase mismatch plays no role 
in the collapse dynamics and we henceforth set (3n = in 
Eq.(|l]). In this case, equations (0) preserve the individual 
powers Pn and, for D >2, the virial integral (0) satisfies 



diK -H 

z - p 



(12) 



where the strict equality still applies to _D = 2 only. From 
Eq. (|l2[) , we infer that collapse inevitably develops when 
H + S^P/4: < 0. This requirement reflects the inter- 
play between the focusing nonlinearities in the invariant 
H and the increase in S'^ of the total mean-square ra- 
dius due to walk-off. As a result, the finite distance 
zf. at which T{z) tends to vanish in case of collapse is 
larger than the collapse focus z° corresponding to S = 0. 
Therefore, the collapse dynamics should be delayed and 



possibly arrested when 6 0. Let us check this prop- 
erty with the (2 -I- l)-dimensional Gaussian beams (|ll|). 
Without walk-off, these beams possess the Hamiltonian 
H = (P/A2)(l - P/Vth) with the collapse threshold 



Vth = 



47r(l + p)2 



1 



2pA' 



(13) 



where p = Pi /P2 . In the presence of walk-off, the esti- 
mate (O) yields the effective critical power for collapse 



Ve =Pth{l 



(14) 



which increases with the walk-off length 5. Figure 2 
shows this critical power for equal waves with Pi = Pj. 
The threshold (^^ is associated with a total collapse 
taking place at the center (r) = and it raises with 
6, which is confirmed numerically. At a given power, 
walk-off is thus able to delay and even arrest such a 
collapse, as it can be guessed from the virial expression 
I{z) = i4zyA^Vth){V! -P) + A\ 
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FIG. 2. Collapse threshold versus 5 = Sx{Sy = 0) for equal 
waves with A=l. The solid curve, where error bars have been 
specified for the sake of accuracy, represents the numerical 
result for P\=P2 and A=2. The dotted lines correspond to 
the analytical predictions Vth{t + 5^), Svr, and 'Pth=4:-K. 



The agreement between the theoretical estimate (|Tj) 
and the numerical simulations is reasonably good for 
small S values. For large 6, Figure 3 displays that walk- 
off significantly competes with the collapse by modify- 
ing the trajectories of both components in the transverse 
plane and by detrapping them. The centroids of the com- 
ponents are then shifted from (f) = and further fuse 
again at the origin, where collapse occurs [Fig. 3(a)]. This 
explains the discrepancies in Fig. 2 between the true col- 
lapse threshold and for a total collapse on the center. 
Alternatively, when the two waves contain enough power 
for promoting individual collapses (Pi = P2 > 47r), walk- 
off can make them separate from each other and collapse 
on their own center of mass, far away from the origin 
[Fig. 3(b)]. In this case, Eq.(p^) emphasizes that walk-off 
becomes important whenever \H\ <C S^P/A, which yields 
a bound for 6, namely S > {2/ A)y/ P/Vth — 1, above 
which the two components can collapse separately with 
individual mean-square radii tending to zero. This re- 
sults in a saturation plateau at P = Stt, towards which 
the collapse threshold saturates from 6 > I. Conversely, 
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in the diffraction regime where collapse does not occur, 
the two components can undergo one or more crossings 
before decaying into noise, or continue to walk off away 
from each other while they both decouple and disperse. 
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FIG. 3. Dynamics of coupled waves undergoing walk-off for 
(a) P=16 witli (5=0.3 (solid), 5=0.7 (dotted), 5=0.8 (dashed), 
5=1.0 (dash-dotted), and (b) P = 40 with S= 1.0 (solid), 
5=1.3 (dotted), 5=1.6 (dashed), 5=3.0 (dash-dotted). 
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To understand the motion of the individual centers of 
(f„(z)) = P-^ J r\un\'^dr, we multiply Eq.(0) by 
and (Vj^u* ) and combine the results to derive 



|u3-n|'Vi|u„prf^^0, (15) 



where dz(r„)(0) = (5„. The integral in (|l^) measures the 
flux induced on the nth wave by its neighbour, which 
can make (r„) periodically displace to a fixed point 
then go back to the center. For equal components mod- 
elled by Gaussians |u„|2 = (P/27rA2)e-l^--<^">l'/^' lo- 



cated symmetrically from the origin, Eq.(]15|) reduces to 
dlirn) + (AP/^A4)(r;)e-2<'-'.>VA= ^ q ^nd (f„) has in- 
deed an harmonic motion, which is damped in collapse 
regimes [A(2) 0] and relaxes otherwise [A(z) +oo]. 
Along one arch of oscillation, we can moreover sup- 
pose that the waves keep a fixed width A, so that the 
same equation provides the relation linking (5„ to the 
maximal displacement of the beam centroid : (rJJ) = 
^^ln(l - 27rA2J2/AP)~i. Hence, the constraint that 

oscillations develop only if (r") remains finite imposes a 
critical value to the walk-off length, namely Sf^ < S'^ — 
AP/2ttA'^. This critical value, below which both waves 
are refiected back to the origin, can also be estimated by 
considering two beams launched in the medium at angles 
6*0 and —6*0 with respect to z, such as Sn = 2no tan 610. By 
following the ray trajectory associated with one beam, ge- 
ometrical optics predicts that this ray experiences a total 
internal reflection if 6*0 satisfies cos^o = '^o/('t-o + '^(0)), 
where n(0) and uq denote the values of the refractive in- 
dex on the axis and at infinity, respectively. This relation 
then expresses as tan^ = ^/"-qj where / is the light in- 



tensity at center. For two beams symmetrically displaced 
from the origin, / is mainly given by the interaction part 
of H, i.e., Hint/P ^ AP/SttA^, which thus restores 6c- 

In conclusion, we have shown that in Kerr media FWM 
lowers the threshold power for collapse of two coupled 
waves, when they are close to exact phase matching 
only. The resulting enhancement of self-focusing, how- 
ever, ceases at moderate mismatch values, for which 
the waves undergo sharp resonances. On the other 
hand, walk-off delays and may stop the self-focusing 
by altering the beam trajectories. For comparison, two 
weakly-separated (l-|-l)-dimensional solitons in birefrin- 
gent fibers are known H to either bounce back and forth 
while forming a trapped state, or separate by developing 
nonlinear oscillations, whenever their amplitude is above 
or below a threshold depending nonlinearly on the walk- 
off length. At higher dimension, analogous behaviors also 
characterize high-power pulses, whose components can 
either mutually attract or separate with walk-off. How- 
ever, the self-focusing dynamics here offers a richer vari- 
ety of interaction regimes, among which, e.g., two cou- 
pled components, being even initially superimposed, can 
split before fusing into a single lobe or reversely collapse 
after splitting. Finally, our analysis suggests that, in sat- 
urable optical media, high-dimensional beams should still 
be capable of coalescing into a robust light bullet, which 
may be desirable for optical switching applications. 
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